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O 
O 

^ I The process of star formation from metal-free gas is investigated by following 

I the evolution of accreting protostars with emphasis on the properties of massive 

pL| I objects. The main aim is to establish the physical processes that determine 

OO I the upper mass limit of the first stars. Although the consensus is that massive 
stars were commonly formed in the first cosmic structures, our calculations show 

^ _ that their actual formation depends sensitively on the mass accretion rate and 

lO ' its time variation. Even in the rather idealized case in which star formation 



is mainly determined by Mace, the characteristic mass scale of the first stars is 



crit — 



m 

CN I rather uncertain. We find that there is a critical mass accretion rate M( 

4 X lO~'^M0yr~-^ that separates solutions with Mace < Merit in which objects 
O I with mass ^ lOOM© can form, provided there is sufficient matter in the parent 

i-^l clouds, from others (Mace > Merit) where the maximum mass limit decreases 

^ I as Maec increases. In the latter case, the protostellar luminosity reaches the 

I Eddington limit before the onset of hydrogen burning at the center via the CN- 

^ I cycle. This phase is followed by a rapid and dramatic expansion of the radius, 

possibly leading to reversal of the accretion flow when the stellar mass is about 
, 100 Mq. 

' Under a realistic time dependent accretion rate that starts at high values (~ 

1O~^M0 yr~^) and decreases rapidly in the high mass regime (M* > 90Mq), the 
evolution follows the case of Maec < Merit and accretion can continue unimpeded 
by radiation forces. Thus, the maximum mass is set by consideration of stellar 
lifetimes rather than by protostellar evolution. In this case, the upper limit can 
be as high as ~ 6OOM0. 

We consider also the sensitivity of the results to the presence of heavy elements 
with abundances in the range Z = 5x 10^^ Zq to 5 x IO^^Zq. The main evolution- 
ary features of protostars are similar to those of metal-free objects, except that 
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the value of Merit increases for metal-enriched protostars. Since the accretion rate 
is lower in a slightly polluted environment, the condition Mace < M^rit is expected 
to be more easily met. We find that for metallicities below ~ 10~'^Zq, where ra- 
diation forces onto dust grains in the flow are negligible, a slightly metal-rich gas 
favors continued accretion and the formation of very massive stars. 

Subject headings: cosmology: theory — early universe — galaxies: formation — 
stars: Population III — stars: formation 

1. Introduction 

The problem of the characteristic mass of the first stars that formed in the universe is 
still largely unsolved. Recent theoretical and numerical studies are revealing in leaps and 
bounds the physical state of star forming regions within the first cosmic structures. In the 
framework of cold dark matter (CDM) cosmology smaller objects tend to collapse earlier. 
In these models, a moderately rare object of ~ 3a overdensity collapses and viriahzes at 
about z ~ 30. However, only sufficiently massive objects can cool in a Hubble time and 
become luminous by forming stars. The minimum cooling mass at 2; ~ 30 is Mmin ~ 10^ Mq, 
and is a decreasing function of the collapse redshift (e.g., Haiman, Thoul, & Loeb 1996; 
Tcgmark ct al. 1997; Fuller & Couchman 2000). After virialization, objects more massive 
than Mmin, called "first objects", can collapse gravitationally owing to molecular hydrogen 
coohng (e.g., Palla, Salpeter, & Stabler 1983, hereafter PSS), and fragment into high-density 
clumps. According to numerical studies, the mass scale of these clumps is rather high, of 
the order of IO^Mq (Abel, Bryan, & Norman 2000, 2002; Bromm, Coppi, & Larson 1999, 
2002; Tsuribe & Inutsuka 2001). However, the formation of fragments of lower mass, down to 
< IMq. is also possible (e.g., Yoshii & Saio 1986; Uchara ct al. 1996; Nakamura & Umemura 
1999, 2001, 2002; Uehara & Inutsuka 2000; Omukai 2001). Indeed, the recent discovery of 
HE0107-5240, a star of 0.8 Mq in the galactic halo with virtually no metals ([Fc/H] =— 5.3), 
attests the ability of very metal-poor gas clouds to give birth to low-mass stars (Christlieb 
et al. 2002). 

After fragmentation, the evolution of gravitationally unstable clumps proceeds in a 
highly no n- homologous fashion, with the central parts collapsing first. This runaway phase 
is induced by cooling provided by H2 line radiation at densities n < lO^'^cm^^, and by 
H2 collision-induced emission at higher densities. The resulting gas temperature is nearly 
constant at several 10^ K, and the innermost region of ~ IM© becomes fully molecular due 
to the three-body reaction (PSS). At densities n ~ lO^^cm"^, the cloud becomes optically 
thick to collision-induced absorption and H2 dissociation works as an effective cooling agent. 
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Finally, without further coohng mechanism, at n ~ lO^^cm"^, a small hydrostatic core of 
mass ~ 10~^ Mq is formed (Omukai & Nishi 1998), similar to the conditions found in studies 
of the collapse of clouds of standard solar composition, i.e., the stellar (or second) core 
according to Larson (1969). 

In the primordial clouds described by the numerical simulations, the outcome of the 
runaway collapse is a tiny (~ 10~^ Mq) protostar, surrounded by a large amount of reservoir 
gas (~ 10'^ ^o)- As a result, the protostar can grow by several orders of magnitude in mass 
by accreting the envelope matter. The mass accretion rate onto the protostar is determined 
by the radial density distribution at the time of core formation, and is related to the prcstcllar 
temperature or, equivalcntly, to the effective sound speed Cg by (e.g.. Stabler, Shu, & Taam 
1980): 

Ma.c ~ |. (1) 

The evolution of zero-metal stars in the main accretion phase was studied by Stabler, 
Palla, & Salpeter (1986; hereafter SPS) up to 10.5 M© and by Omukai & Palla (2001; here- 
after. Paper I) for higher masses under a constant mass accretion rate of 4.41 x IO^'^Mq yr~^, 
corresponding to a gas temperature of 1700 K in eq. (1). 

The mass of the forming star is set at the time when protostellar accretion stops. In 
paper I, we have found that fast accreting protostars enter a phase of rapid expansion at a 
mass of ~ 300 Mq when the luminosity becomes very close to the Eddington limit. This 
event may determine the onset of a powerful stellar wind driven by radiation pressure that 
can effectively quench further accretion of circumstellar material. In that event, after a 
transient internal readjustement, the protostar is expected to settle as an ordinary zero-age 
main sequence (ZAMS) star. 

The upper limit of ~ 300 Mq to the protostellar mass has been obtained using a con- 
stant value of the mass accretion rate. In realistic collapse calculations, deviations from the 
isothermal assumption are expected to occur, and even in the ideal case the mass accretion 
rate is found to vary in time (e.g., Foster & Chevalier 1993; Tomisaka 1996). Therefore, it 
is important to extend our initial study to a larger set of conditions in order to explore the 
robustness of the results described above. Along with differences in the mass accretion rate, 
we will also address the issue of the impact of a small abundance of heavy elements. It is 
expected that metals may alter the mass hmit drastically as a result of the existence of a 
dust shell and because of steller evolution effects. 

In §2, we sketch our numerical approch. Then, in §3 we describe the results of the 
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main accretion phase obtained with a constant and time variable Mace- The effects of a 
finite, non-zero metallicity on protostellar evolution are discussed in §4. In §5 we discuss the 
implications of our work. Finally, in §6, we give the conclusions. 

2. Numerical Approach 

The basic strategy and equations used here are the same as in SPS. In that scheme, 
the protostellar evolution is treated as a sequence of steady state accretion flows onto a 
growing core (sec Figure 1 of SPS for a definition of the various regions). The core is 
assumed to be in hydrostatic equilibrium and the ordinary stellar structure equations are 
applied. If the gas in front of the accretion shock is optically thin, no envelope model is 
constructed and the boundary condition given by cq. (4a) of SPS is applied at the core 
surface. In case the preshock gas is optically thick, we integrate the equations inside the 
radiative precursor from the photosphere to the core surface (see eqs. 7a-d of SPS). Outside 
the photosphere, we assume a free-falling flow. The main difference with the SPS work 
is in the use of updated opacity tables: for T < 7000 K, we take the results of Lenzuni, 
Chernoff, & Salpeter (1991) for a gas composition X = 0.72, y = 0.28, and the OPAL 
opacity at higher temperatures (Iglesias & Rogers 1996), with a slightly different hydrogen 
abundance X = 0.70). For partially polluted gas, we have used the results of Alexander & 
Ferguson (1994) with X = 0.70, Y = 0.30 for T < 7000 K, and the OPAL opacity for higher 
temperatures. Although the primordial He abundance is about Y = 0.24 (e.g., Galli et al. 
1995; Izotov & Thuan 1998), we used rather high value of F > 0.28 as stated above because 
of the lack of opacity tables for primordial gas. This causes no significant error. 

The calculations are started at the protostellar mass of O.OIMq. The initial models are 
arbitrary and constructed using an entropy distribution of the core given by 

s{M) = So + /3(M/O.O1M0)2 , (2) 

where M is the enclosed mass within a spherical STirfacc, Sq = —13.05 and f3 = 7.488 in 
units of TZ, and the zero-point is arbitrarily set at T = 2.05 x 10^ K and p = 5.16 g cm~^. 
With the mass accretion rate used by SPS, M^cc = 4.41 x 10~^ Mq yr~^, the initial central 
temperature is 10^ K. This value increase (decreases) for higher (lower) accretion rates. The 
initial conditions do not affect the later evolution, since the protostar adjusts quickly during 
the "decay-of-transient phase" (for the protostellar mass < O.IM©; SPS) to the state 
appropriate to an accreting object. 



3. Formation of Primordial Stars: 
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Effects of Variation of the Mass Accretion Rate 

3.1. Evolution under different mass accretion rates 

First, we study the effects of using different, but constant values of Mace on the evolution 
of primordial protostars. We define Mfid = 4.41 x 10~^ Mq yr~^ as the fiducial value of the 
accretion rate used by SPS and in Paper I. To follow the dependence of the evolution on 
-^acc, wc have studied four cases, varying Mace from 1/4 Mad to 2Mfid, thus covering the 
range 1.1 - 8.8 x 10"^ M© yr'^ 

The evolution of the protostellar radius as a function of the accretion rate is shown 
in Figure 1. As the calculations start from rather arbitrary initial conditions, the protostar 
adjusts itself to the state appropriate to an accreting object by entropy redistribution. As 
a result, the fast entropy variations cause the bumps in radius seen in Figure 1 at masses 
<0.03 Mq. After this minor transient phase, the evolution proceeds smoothly and the 
similarity of the curves allows the identification of three distinct phases. The first two are 
characterized by a common trend where the protostellar radius i?* increases almost linearly 
with the protostellar mass M* in the logarithmic scale, reaches a peak between ~10 and 
20 Mq, and then drops substantially. The adiabatic accretion phase (Phase I, M* < a 
few to ~ IOM0) is marked by a gradual expansion of the radius and lasts until the arrival 
of the internal luminosity wave at the surface, whose eruption causes the sudden swell of 
the external layers (e.g., Palla & Stabler 1991). In the Kelvin-Helmholtz contraction phase 
(Phase 11, a few-lOA'/Q < M* ^ 30 — 60Mq) the increased gravitational pull of the growing 
star stops the expansion and forces the star to contract toward the conditions appropriate 
for a main sequence object. The third phase varies markedly, depending critically on the 
adopted value of Mace- 

If the rate is comparable to Mgd radiation pressure causes the abrupt expansion of the 
surface layers (Phase 111b) seen in the two upper curves of Figure 1 computed for Mace = 
(1 — 2) Mfid- The protostellar radius increases tremendously and the impulsive release of 
energy likely causes the stripping of the external layers of the protostar by a radiation- 
driven stellar wind. As shown in Paper I, this event signals the end of the accretion phase 
and sets the maximum value of the protostellar mass. In the opposite case where Mace < -^fid) 
contraction proceeds unimpeded by accretion. In the cases shown in Figure 1 and computed 
for Mace — 1/2 and l/4Mfid, the protostar reaches the ZAMS at a mass of ~ 50 Mq (Phase 
Ilia). 

The global evolution of the core interior is shown in Figure 2 for the low mass accretion 
rate (Mace = l/4^fid) and in Figure 3 for the case Mace = Mgd- In these figures, the bottom 
panels illustrate the evolution of the central variables, as well as those of the position of the 
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temperature maximum for the same cases. 

The overall evolution of the interior luminosity is shown in the top panel for Mace — 
l/4Mfi(j and in the bottom panel of Figure 4 for Mace — 2Mfid- The contribution to the total 
luminosity provided by nuclear reactions of deuterium and hydrogen (through the pp-chain 
and CN-cycle) at the center of the protostar is also indicated. Notice how the onset of each 
burnig stage is shifted to higher masses for bigger values of Mace- 

Having sketched the main features of the evolution, we now describe in detail the results 
of each evolutionary phase identified above. 



3.1.1. Phase I: adiabatic accretion and propagation of the luminosity wave 

Early in the evolution, the opacity in the interior is due to free-free absorption (/tg cx 
pT~^-^) and is very high due to the low temperatures (T < 10^ K). As shown in Figure 4, 
the luminosity is extremely low and consequently the cooling time of the internal regions 
(i.e., the Kelvin- Helmholtz time = GM^/R^L^) is longer than the evolutionary timescale 
(i.e., the accretion time tace = M*/Macc)- As a consequence, the accreted material piles up 
without further cooling, after passing through the accretion shock and the settling region. 
During this adiabatic accretion phase, the protostar swells gradually as the mass increases. 
The protostellar radius obeys the relation (see equation 21 of SPS) 

i?.ocMr^M^^\ (3) 

that can be derived assuming that the luminosity is dominated by the accretion contribution 
and that the opacity in the envelope is due to H~ bound- free absorption, scaling as k oc T^^-^. 

The increase of the temperature accompanied by the growth of the protostellar mass 
produces a significant decrease of the opacity. Eventually, the heat contained inside the 
protostar begins to propagate outward as a luminosity wave (see Figure 10 of SPS). Then, 
the interior luminosity increases dramatically at a mass of ~ 4M0 for Mace = l/4Mfid, 
and at ~ SMq for Mace = Mid (see Figure 4). The arrival of the luminosity wave at the 
surface produces the sudden sweUing of the protostar, marking the end of Phase I. In the 
course of this expansion, a short-lived surface convective region develops which extends to 
about 20% in mass (Figs. 2, 3). 

In the case of low Mace (Figure 2), we can see a temporary, but active phase of deuterium 
burning around 3M0, when the nuclear energy generation rate exceeds 0.1L*/M*. This 
episode is just the consequence of the low value of L* before the arrival of the luminosity 
wave (see the top panel of Figure 4). However, the deuterium luminosity is very small 
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and deuterium is hardly consumed with no significant consequences on the evolution. 

As a result of the outward increase of the entropy, the temperature is highest off-center, 
at an internal mass of M/M^ ~ 0.1. The position of the off-center temperature peak is shown 
in the top panels of Figs. 2 and 3 and is maintained throughout the adiabatic accretion phase. 



3.1.2. Phase II: Kelvin-Helmholtz contraction 

After the relaxation of the interior and the passage of the luminosity wave, the major 
opacity source in protostars is provided by electron scattering, except in a small central 
region where free-free absorption is dominant, and in a thin surface layer characterized by 
bound-free absorption of H and He. 

The condition of radiative equilibrium and of an opacity independent of density and 
temperature throughout most of the interior determines the well known mass-luminosity 
relation L* oc M^, and the Kelvin-Helmholtz time becomes 

t - oc ^ (4) 

If the radius is sufficiently large, txH < ^acc and the core shrinks due to fast energy loss, 
while in the opposite case the core swells adiabatically. The radius of the protostar actually 
adjusts towards equality of the two timescales txH — ^acc, leading to the relation 

oc M^c/Ml (5) 

Small deviations from this relation are caused by minor contributions to the opacity from 
processes other than electron scattering. 

Figure 5 shows the evolution of the interior luminosity to mass ratio, L^/M^,, as a 
function of the protostellar mass M^,. During KH contraction, L^/M^ increases rapidly and 
the mass-luminosity relation follows approximately a unique relation independent of the 
value of the mass accretion rate, as expected in the ideal case of constant opacity. Also, the 
photospheric luminosity, which is the sum of the interior (L*) and accretion (Lace) luminosity, 
is independent of Mace; since Lace oc Maee/-R* a-nd R^ oc Mace- 

At the beginning of contraction, the maximum temperature exceeds lO^K, and deu- 
terium can start burning. Because of the off-center temperature peak, burning propagates 
both inward and outward, scouring out all of the available deuterium (see Fig.s 2 and 3). 
Under the high accretion rates considered here, however, D-burning plays a minor role in 
the evolution and contributes at most ~l/4 of the interior luminosity (see Figure 4). The 
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off-center temperature peak is only a transient property of the core and disappears at 20Mq 
for Mace = l/4Mfid, and AOMq for Mace = -^fid, as seen in Figs. 2 and 3. 

After deuterium is completely destroyed and the protostar has contracted further, the 
pp-chain becomes operative in the center. However, even in this case, the energy generated 
by these reactions is not enough to contrast contraction effectively, owing to the saturation 
of the energy generation rate at high temperatures (see Figure 4) . Therefore, despite nuclear 
heating, protostars can continue the phase of gravitational contraction. 

3.1.3. Phase Ilia: Settling onto the ZAMS with continuing accretion (M^cc < l/2Mfidy) 

For accretion rates of Mace = 1/2 and l/4Mfid, during contraction the central tempera- 
ture becomes high enough to synthesize carbon via the CN-cycle and H-burning becomes a 
significant source of luminosity. Figure 5 shows that at this time gravitational contraction 
is halted and the interior luminosity to mass ratio L^/M^ stops increasing. The protostar 
relaxes quickly to the ZAMS and the models converge toward a single mass-luminosity-radius 
relation, as shown in Fig.s 1 and 5. With H-burning working as a thermostat, the central 
temperature remains almost constant at about 10^ K (see the bottom panel of Figure 2). 

At this point, the interior luminosity is essentially supplied by nuclear reactions (see 
the top panel of Figure 4). Since part of this energy source is used to heat the interior, 
the specific entropy in the convective core increases as the stellar mass grows. This results 
in the decrease of the central density, while the temperature remains constant. The mass 
fraction of the convective core increases (see the top panel of Figure 2) and the CN abundance 
reaches very quickly a level of ~ 10~^, thereby gradually increasing (see the bottom panel 
of Figure 2). 

As shown in Figure 5, throughout the evolution the luminosity remains below the Ed- 
dington value by a rather large margin. The latter has been computed using electron scat- 
tering as the main source of opacity in the accretion flow. As a result, accretion can continue 
without experiencing the effects of the retardation force exerted onto the ionized radiative 
precursor. Therefore, the upper limit of the protostellar mass formed in these conditions is 
set by the amount of circumstellar material left during accretion and not by the feedback 
effect of radiation pressure. 
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3.1.4- Phase Illb: Rapid expansion due to radiation pressure (Mace > -^fidj 

For high accretion rates, the protostar undergoes a phase of violent expansion as soon as 
the interior luminosity becomes comparable to the Eddington limit. This process is caused 
by the strong radiation force exerted both on the stellar surface and the radiative precursor. 
The result is that fast accreting protostars never converge to an ordinary ZAMS structure. 
Let us see more in detail how this condition is achieved. 

During the contraction phase, the luminosity-to-mass ratio increases until the condi- 
tions at the center are appropriate to start nuclear burning (see Fig. 5). This happens at 
progressively large masses as the accretion rate increases. For sufficiently high Mace, the 
total luminosity, which is the sum of the interior and accretion luminosity, reaches the Ed- 
dington hmit before the ignition of hydrogen. This is the case for Mace — 2Mfid, while for 
the lower value of Mace = Mfid the situation is more complicated. In the accreting envelope, 
gravitational acceleration is approximately balanced by the deceleration due to the radia- 
tion force. Then, the ram pressure onto the protostellar surface decreases suddenly, and the 
protostar begins to expand without violating the Eddington limit by lowering the accretion 
luminosity Lace oc R^^. During expansion, the temperature drops, while the radiation force 
and the opacity increase at the surface, the latter peaking sharply owing to the contribution 
of ionization of trace amount of atoms. This is mostly a surface effect since the main opacity 
source is still electron scattering both inside the protostar and in the accreting envelope. 
By this mechanism, the expansion is accelerated and possibly results in the stripping of 
the surface layers and in the reversal of the material in the accreting envelope. Thus, the 
main phase of accretion is over and the protostar can rapidly relax onto a ZAMS structure. 
According to the results shown in Fig. 1, this critical phase is achieved at a mass of ~ 90Mq 
for Mace = 2Mfid, and ~ SOOMq for Maec = Mgd. Therefore, our models show that the upper 
limit of the mass shifts toward smaller values for progressively higher values of the mass 
accretion rate. 



3.1.5. A general consideration 

We have seen that, depending on the actual value of the accretion rate, the final fate 
of primordial protostars can be quite different and leads to strongly different values of the 
maximum limit of their mass. This suggests to look deeper into the physical origin of the 
bifurcation displayed in the mass-radius relation of Figure 1. 

The total luminosity of a protostar, Ltot, consists of the interior, L^,, and accretion 
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luminosity, Lace: 

Ltot ^L^ + ivacc ^L^ + GM^M^/R^. (6) 

Imagine that the protostar reaches the ZAMS star while continuing accretion, as in the cases 
of Mace < l/2Mfid- Since a unique mass- luminosity-radius relation holds for ZAMS stars, 
the total luminosity would be 



^tot 



^ ivzAMS + GM.M^^/RzAMS , (7) 



where we have set i?* = -Rzams(-M*) and = -^zams(-M*). If the accretion rate is larger 
than the critical value 

47rcil!zAMS -^ZAMS^ ^ _ in-3 



Mc,H = ^ 1 - ^ 4 X 10-^Me/yr , (8) 

where Kes is the electron scattering opacity, the total luminosity exceeds the Eddington 
limit ivEdd, implying that the star cannot adjust to the ZAMS structure at these accretion 
rates. Although the critical accretion rate is a function of the mass of the central star, the 
dependence is weak considering the — — relations. 

Note that the critical accretion rate Merit is the value where the total luminosity L^t 
reaches L^dd at the same time as the protostar reaches the ZAMS. Thus, we can derive 
Mciit by using the ZAMS quantities, although protostars start expanding before reaching 
the ZAMS for Mace > Merit- For Mace < Merit, -^tot remains below Lsdd up to the point 
where the protostar reaches the ZAMS, while for Mace > Merit, -^tot equals L^dd in the 
course of the Kelvin-Helmholtz contraction and the protostar starts expanding. 

The weak dependence of Merit on the protostellar mass M* can be seen more clearly if 
we define the Eddington radius as the radius of the core where the total luminosity (defined 
in eq. 7) equals the Eddington limit L^dd- 

p _ ^esMacc /qx 

- Anc{l - r) ' 

where F = L^/Lsdd- 

In the two panels of Figure 6, we show the Eddington radius and the core radius for 
the four values of the accretion rate. In all cases, when the core mass is small and F ^ 1, 
the Eddington radius is constant. Later on, the evolution differs. Consider first the case 
for low Mace- During the contraction phase, the luminosity-to-mass ratio increases, the core 
shrinks, but the Eddington radius starts a slow rise. However, i?Edd never climbs up to i?*, 
since the core meets the ZAMS conditions before and its radius also increases as a result of 
the heating provided by nuclear burning. Thereafter, i?Edd remains below the ZAMS radius 
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and accretion can proceed unimpeded by radiation forces. Since the i?Edd/-R* ratio remains 
approximately constant, the critical accretion rate is almost independent of the protostellar 
mass. 

If, on the other hand, the protostar reached the ZAMS structure with an accretion flow 
Mace > -^crit, the total luminosity would exceed the Eddington limit. Unlike the previous 
case, at some point of the evolution the Eddington and core radii would become equal (see 
the bottom panel of Fig. 6). The core now undergoes the runaway expansion, keeping the 
luminosity slightly below the Eddington limit. The increased radius reduces the accretion 
luminosity, thus allowing more gas to flow for a short period. 

3.1.6. Properties of the radiative precursor 

Except for the short-lived phase of the peak radius around IOA/q, primordial protostars 
are embedded within an optically thick envelope. The spatial extent of the precursor, which 
is the optically thick part of the accretion flow, corresponds approximately to that of the 
ionized region. However, the ionized region around the accreting protostar differs from an 
ordinary HII region in that it is opaque to optical photons, and both matter and radiation 
are in the thermodynamic equilibrium. 

In Figure 7, we show the evolution of the photospheric radius for two cases with Mace = 
l/4Mfid and Mgd. During the adiabatic phase, most of the luminosity is generated at the 
accretion shock and the ratio of the photospheric to core radius remains nearly constant, 
with a value of ~1.4 (see also eq. 19 of SPS). After the luminosity wave reaches the surface, 
the interior luminosity exceeds the accretion contribution. As a result of core contraction 
and of the increased heat deposit in the precursor, the ratio of photospehric to core radius 
increases dramatically. The trend is reversed during the final stages of the evolution. 

The evolutionary tracks of the photosphere in the HR diagram are displayed in Fig- 
ure 8 for the four accretion histories. For comparison, we also show the evolution of the 
protostellar surface and locus of the ZAMS for metal-free stars. The distinction between the 
surface temperature and the ZAMS effective temperature after reaching the MS phase is a 
consequence of the difference in boundary conditions at the stellar surface. Because of the 
optically thick envelope, the hot surface of the protostar remains invisible in the optical. The 
photospheric temperature is locked at ~6000 K due to strong sensitivity of the H~ bound- 
free opacity on temperature. However, for low Mace, the tracks make an abrupt shift toward 
higher temperatures for M* > lOOM©. As the protostellar mass increases, the density in 
the accreting envelope decreases, while the temperature and ionization degree both increase. 
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Electron scattering replaces the H" b-f absorption as the major source of opacity at the pho- 
tosphere. Since this process is independent of temperature, the effective temperature rises 
as the density in the accreting envelope decreases. Correspondingly, the photospheric radius 
decreases sharply (Figure 7). The tracks would then join the corresponding ZAMS models 
at masses exceeding lOOOM©, but we cannot determine exactly where this would occur since 
such models are not available in the literature. 



3.2. Protostellar Evolution with a mass dependent accretion rate 

The results described so far have been obtained under the rather limiting assumptions 
that the accretion rate is constant throughout the evolution and of large magnitude (> 
10~^MQyr~-'^), at least two orders of magnitude more than what is normally considered in 
the coUpase phase of present-day molecular cores. That Mace is so large depends on the 
thermal and chemical properties of the primordial gas, which is known to be a very poor 
radiator owing to the lack of molecules and dust grains. Since the gas temperature during 
collapse is mainly determined by the excitation of H2 molecules, it is natural to expect that 
collapse models find typical values well in excess of several hundred degree, thus justifying 
the large values of Mace through eq. (1). On the other hand, the use of a single value of 
the accretion rate is less motivated since, as we have noted in §1, numerical simulations 
have shown a rather strong dependence on time of the amount of material that can actually 
collapse, even in the case of strict isothermality. 

Thus, it is natural to ask how robust is the sequence of events that we have described 
in the previous sections in view of a more general behavior of the gas dynamics, and hence 
of the time history of the accretion rate. Since the key factor that determines the ultimate 
fate of massive protostars is whether Maec is above or below the critical value Merit — 
4 X 10~^ Mq yr^^, we should consider how and when this situation is met under more 
realistic conditions. 

A useful guide is provided by the recent 3D hydrodynamical simulations of the formation 
of the first stars presented by Abel, Bryan & Norman (2002; hereafter ABN). Starting from 
cosmological initial conditions, these authors have followed the evolution of a collapsing 
primordial clump up to the stage when optical depth effects become very large, i.e., when the 
central number density reaches ~ lO^^cm"^. Due to numerical limitations and their optically 
thin approximation for H2 line radiation, ABN could not follow correctly the subsequent 
history of protostar formation and accretion phases. However, they provide the run of the 
accretion time as a function of enclosed gas mass which can be converted into an effective 
accretion rate that depends on the growing protostellar mass. Using the results shown in 
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Fig. 5 of ABN, we obtain the following fit 

r 7.76 X 10-3 (M,/Mq)-o-24 Mq yr-i 0.6 Mq < M, < 13 Mq 
Mabn = < 1.51 X 10-3 (M,/M0)O-^i Mq yr-^ 13 < M, < 60 (10) 
i 11.2 (M^/Mq)-!-^^ yj.-i 60 Mq < M,. 

In order to avoid computational difficulties, we have used Mace = 8.8 x lO-^M© yr-^ for 
lowest mass models (M* < 0.6 Mq). Thus, the accretion rate starts at very high values, 
it remains well above ~10-3 Mq yx~^ up to about 60 Mq, and then falls rapidly for very 
massive stars. Most importantly, Mabn drops below the critical value for M* > 95 Mq. 

The evolution of the protostar with the accretion rate given by equation (10) is displayed 
in Figure 9. For comparison, we also show the results at fixed Mace discussed before. In 
Figure 10, the evolution of the photospheric radius is compared to that of the core. During 
the initial adiabatic phase, the protostar expands more slowly than in the constant Mace 
case in response to the decrease of Mabn for M* < 13 Mq. The thermal relaxation and the 
appearance of the luminosity wave with the expansion of the radius occur at approximately 
the same mass as for Mace = ^fid) since the instantaneous value of Mabn is nearly the same 
as the fiducial one (0.97 Mgd at 13 Mq). 

During Phase II, the protostar contracts more gradually than in the fiducial case since 
Mabn reaches a minimum at 13 Mq and increases thereafter. The core shrinks down to 
its minimum value of ~ 20 Rq at about 50 Mq and then rebounds when radiation pressure 
effects begin to appear. This is the most delicate part of the evolution that depends critically 
on the behavior of Mabn- As Fig. 9 shows, the phase of violent expansion is avoided in this 
case and the protostar resumes further contraction. At about 90 Mq, hydrogen is ignited 
at the center and a convective core begins to develop. Then, the accretion rate falls below 
the critical value at 95 Mq. Shortly after, the protostar relaxes to a ZAMS star, without 
suffering from surface expansion or a blow out of the accretion envelope. Therefore, accretion 
is expected to continue unimpeded to very high masses. 

With continuing accretion, the final mass of the forming star is set by how much 
material can be accreted during the stellar lifetime (e.g., Larson & Starrficld 1971). Be- 
cause of the sharp decrease of the accretion rate, the accretion time tacc = M*/Mabn — 
3 X 10^(M,/100MQ)2-76yr (for M, > 60Mq) exceeds the steUar fifetimc (~ 3 x 10^ yr) at 
M* ~ 500Mq. Therefore, even if accretion goes on throughout the entire stellar lifetime, 
the stellar mass would not exceed 500Mq by a large margin. A more detailed analysis has 
concluded that if the accretion continues, the mass attained during the entire stellar lifetime 
is about 600Mq (ABN). In reality, when the accretion rate drops below a certain value, the 
star is hot enough to produce a copious amount of ionizing photons that create an ordinary 
HII region, whose expansion can effectively reverse the infalling gas. Omukai & Inutsuka 
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(2002) have shown that the mass hmit set by this process is very similar to the hmit imposed 
by the main-sequence hfetime (< 46OM0; see their eq. 27). 

Returning to the evolution of the photospheric radius, in the high mass regime 
6OM0), the rapidly decreasing accretion rate causes the abrupt turn around shown in Fig. 10, 
after reaching the peak value of ~ 1000 it!© at about 100 Mq. At a mass of ~ 500 M©, the 
radiative precursor disappears and the accretion envelope becomes optically thin: At this 
point the hot stellar surface becomes directly visible from outside. 

4. The formation of second-generation stars: 
Effects of metal enrichment 

Second-generation stars are formed out of material that has been contaminated by 
metals released from Population III SNe. Heavy elements affect the chemistry and dynamics 
of the gas in several ways. For example, detailed models have shown that a metal abundance 
of ~ 10^^ Zq can substantially change the cooling properties of gas, thereby lowering the 
minimum mass of cosmological objects (e.g.. Fall & Rees 1985). A metal abundance as low 
as ~ 10^'^ Zq can influence the fragmentation mass scale of star-forming clouds (Omukai 
2000; Bromm et al. 2001a; Schneider et al. 2002). Protostellar evolution is also altered 
through important changes in the opacity and in the initial abundance of the CN elements. 
Therefore, it is important to consider how even a minute enrichment of heavy metals modifies 
the picture of star formation that we have been describing so far. 

Since the preceding discussion has demonstrated that the qualitative evolution of ac- 
creting protostars is the same for fixed and variable mass accretion rates, we have run several 
models at constant Mace for nonzero metallicities. Figure 11 shows the resulting mass-radius 
relation obtained with M^^ = 1/2Mm = 2.2 x 10"^ Mq yr"^ for Z = 10"^ = 5 x 10'^ Zq 
and 10"*^ = 5 X along with the standard zero-metal case. The evolution is basically 

the same, except for the temperature at which the CN-cycle starts operating. At higher Z, 
significant burning begins at lower temperatures. Therefore, the contraction phase is less 
pronounced, yielding a larger core radius on the main sequence. This explains the shift of 
the minimum radius in Figure 11 at higher Z. On the other hand, the protostellar lumi- 
nosity at a given mass is the same in all models. The main effect of metals consists of an 
increase of the critical accretion rate that varies from Merit = 4 x 10"^ Mq yr~^ for Z = to 
9 X 10~^ Mq yr"^ for Z = 10~^ = 5 x 10"^^©. The change of a factor -2 in Merit reflects the 
similar variation of the protostellar radius (through equation 8). Of course, a higher value 
of the critical accretion rate implies that protostars can keep accreting without entering the 
phase of violent expansion. 
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The effect of metals is much more dramatic if we consider other values of Mace- Figure 12 
displays the mass-radius relation computed with Mace = Mfid = 4.4 x 10~^ Mq yr~^ for the 
three metal compositions. The accretion rate is now slightly higher than the critical value 
for the mctal-frcc case, but smaller than that aX Z = 10^^ = 5 x IO^^Zq (6 x 10^^ Mq yr~^). 
As a result, the violent expansion that occurs at Z = is no longer found, even in the lowest 
metallicity case. 

Thus, we conclude that the presence of metals, by increasing the critical accretion rate, 
favors the formation of very massive objects. This result is not intuitive since metals, locked 
up in grains, arc usually believed to reduce the upper mass limit because of the enhanced 
radiation force onto the dust shell surrounding the protostar (e.g., Kahn 1974; Wolfire & 
Cassinelli 1987). To understand this behavior, let us remind that the mechanism that should 
prevent infall is the violent expansion of the surface layers caused by the radiation force on 
the ionized radiative precursor. In the case of high Z, radiation pressure is imparted onto 
the dust shell which is located at a much greater distance from the precursor. 

We can derive a simple estimate of the critical metallicity that separates the two regimes. 
The radiation force on the dust shell becomes bigger than that on the ionized precursor if 
the opacity in the shell Kdust = '^dust('^/'^0)' where K^^g^ ~ 30 cm^ g"-*^ is the present value 
of dust opacity (e.g.. Beech & Mitalas 1994), exceeds that due to electron scattering in the 
ionized precursor, n^rcc = 0.4 cm^ g^^. Here, we have assumed that the dust properties (e.g., 
composition, size distribution, etc.) do not change with metallicity. In this case, the critical 
value of metallicity is given by 

/ K^^^ \ 

Z>Zdr = 0.01Ze ''f _^ 

\ 30cm^g ^ / 

Note that Wolfire & Cassinelli (1987) argued that the dust must be depleted by a factor of 
~ 4 from the average interstellar value for the formation of massive stars. If so, the opacity 
will be lower by the same factor, and this will affect the critical metallicity Zdr as well. 
For Z > Zdr, the infall is reversed by blowing out the outer envelope before the luminosity 
reaches the Eddington limit due to electron scattering. In this regime, metals act against 
prolonged accretion. For Z < Z^t, the radiation force onto the outer envelope is always 
below gravity, since the protostellar luminosity never exceeds the Eddington limit. When 
~ -^^Edd, the radiation force on the radiative precursor and on the surface layers causes 
the violent expansion of the protostar. In this regime, metals work in favor of continuing 
accretion by increasing the critical accretion rate. 
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5. Discussion 

Much of the recent attention on primordial stars has focussed on the properties of mas- 
sive objects. Several authors have argued that the first generation of cosmic structures 
formed mainly massive and very massive stars, and indirect evidence for a top-heavy IMF 
includes chemical abundances and supernova rates in our Galaxy, as well as star formation 
rates in high redshift galaxies (e.g., Larson 1998). Also, the apparent lack or dearth of stars 
in the galactic halo with metallicity less than ZjZ^ ~ 10~^, coupled with the theoretical 
prediction of a critical metallicity for efficient fragmentation at essentially the same vahie 
of Z (Bromm et al. 2001a) have led to the suggestion that the IMF at early times differed 
substantially from its present form. Finally, the initial chemical enrichment of our Galaxy 
and the abundance pattern in the intracluster medium have been considered as compelling 
evidence for a generation of stars with mass >100 M© (Schneider et al. 2002). In particular, 
these authors have postulated that in order to explain the apparent average minimum metal- 
licity of ZjZo, ~ 10~^, the first stars should have had masses in the rather narrow interval 
between 140 Mq and 260 Mq, ending their hves as pair- unstable supernovae and polluting 
the pristine gas with the right amount of heavy elements. Higher and lower mass stars would 
in fact produce massive black holes, without releasing metals (e.g., Heger & Woosley 2002). 

While there are good theoretical reasons to believe that massive stars were common in 
the early universe, the question of their actual formation has not been studied in detail. As 
we know from studies of the local universe, forming massive stars is a rare event that requires 
special conditions and it is not clear that inefficient fragmentation, lack of metals and fast 
mass accretion are per se sufficient ingredients to justify the expectations. The discovery 
of an extremely metal-poor halo star of ~solar mass with an iron abundance [Fe/H]=— 5.3, 
twenty times lower than the previous limit of Z/Zq fa 10~"^ (Christlieb et al. 2002), clearly 
indicates that low-mass star formation in such extreme environments is more pervasive than 
otherwise believed. 

The results of our calculations have shown that the final outcome of the accretion process 
is a sensitive function of the mass accretion rate and of its time variation. Notwithstanding 
the role of other important physical processes that we have neglected (presence of circumstel- 
lar disks, turbulence and magnetic fields), even in the idealized case in which star formation 
is mainly determined by Mace leads to uncertain predictions on the characteristic mass scale 
of the first stars. The existence of a critical mass accretion rate Merit — 4 x 10~^ Mq yr~^ 
separates solutions in which very massive objects with mass ^ 100 Mq can form, provided 
there is sufficient material in the parent clumps, from others where the upper mass limit de- 
creases as Mace increases. In the case of the time dependent accretion rate from the numerical 
simulations of Abel et al. (2002), protostellar growth can continue unimpeded by radiation 
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forces onto the radiative precursor. However, slight deviations of the actual accretion rate 
from the empirical form of eq. (10) may result in a large difference in the final mass of the 
formed stars. Obviously, more refined models of the collapse phase should be carried out to 
verify the robustness or uniqueness of the behavior of the accretion rate in the high-mass 
regime (M* > 60 Mq) which is critical for prolonged accretion. If the decrease of Mace is 
slower than that given by eq. (10), accretion can be stopped at around a few tens to 100 
Mq, thus not allowing the formation of the objects which are considered responsible for the 
initial enrichment of the primordial gas. 

Another important result of our study is the appreciation of the sensitivity of the evolu- 
tion on the metal abundance. In the Z = case, the value of the actual accretion rate from 
eq. (1) is similar to the critical one. For the metal enriched case, Ma.cc is expected to be 
smaller, as a result of the reduced gas temperatures of the prestellar clumps. On the other 
hand. Merit increases at higher Z, as shown in §4. Hence, in a slightly metal-enriched gas. 
Mace would be lower than Merit, leading to continued accretion and favoring the formation 
of very massive stars. 

Considering the mass scale for cloud fragmentation, recent studies have found a tran- 
sition from massive (~ 10^ Mq) to low-mass (~ IMq) clumps at a metaUicity of Z^ag ~ 
10~'^Zq (Bromm et al. 2001a; Schneider et al. 2002). However, as discussed in the previous 
section, there appears to be another critical metallicity. U Z > Z^r ~ O.OIZq, radiation pres- 
sure onto the dust shell surrounding a protostar becomes important and limits the amount 
of matter that can be accreted. Therefore, one can argue that there arc three regimes that 
determine the maximum value of the stellar mass: for Z < Zf^^g-, the limit is determined by 
the the stellar lifetime and/or by the expansion of the ionized region; for Z^ag < Z < Z^^ by 
the fragmentation scale, and by the reversal of the gas inflow due to radiation pressure onto 
the dust shell for Z^j. < Z. In the first case, massive stars form predominantly, while in the 
other two cases their production is possible, although not typical. 

In this work, we have found that if Mace > Merit, radiation pressure on the ionized 
outer layers of the protostar causes a violent expansion of the radius. Since we have treated 
protostellar evolution as a series of hydrostatic configurations, the final outcome of this 
dynamic phase could not be fully investigated. We have suggested that the rapid swelling 
leads to the ejection of the external layers, to the onset of powerful mass loss through a 
wind, and eventually to the complete halt of accretion. Note that stellar pulsations due to 
the e-mechanism triggered by the onset of the CN-cycle can lead to moderate mass loss in 
metal-rich stars, but in the case of metal-free stars the loss rate is strongly reduced (Baraffe, 
Heger, & Woosley 2001). 

There are, however, other possibilities that need to be considered. For example, the 
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stellar wind may not be strong enough to blow away all the envelope material. With no 

accretion, the central star would quickly contract to the main sequence, decreasing sub- 
stantially its mass loss activity. Then, the remaining circumstellar matter would be able to 
collapse and recurrent accretion might be possible. 

Deviations from spherical symmetry is also crucial. The stellar wind might blow in a 
bipolar jet-like configuration, without affecting accretion from the disk. Even if the accretion 
rate through the disk is higher than the critical rate, the protostar would adjust itself in a 
state marginally below Afcrit by shedding the remainder material as a jet-like wind. Such 
energetic jet could disrupt the entire parent clump and quench the supply of matter to the 
disk (e.g., Nakano, Hasegawa, & Norman 1995). 

When the accretion proceeds in a disk-like fashion, the star loses energy faster than 
in the case of spherical accretion since the stellar radiation can escape in the optically thin 
polar directions. The reduced entropy level in the star results in a smaller protostellar radius 
(Palla & Stabler 1993; Hartmann, Cassen, & Kenyon 1997). Thus, at the same protostellar 
mass, the accretion luminosity is higher in the disk accretion than in the spherical case. 
This suggests that accretion in the polar direction is halted earher on. On the other hand, 
if the parent clump is not disrupted by sweeping-up of the matter in the polar direction, 
accretion can proceed through the disk since only the matter at large angles from the disk is 
subject to the radiation force (Nakano 1989). To decide whether the disk-like configuration 
is more suitable or not for continuing accretion, 2D or 3D hydrodynamical studies must be 
carried out, similar to those now available for the present-day case (e.g., Yorkc & Sonnhalter 
2002). Recently Tan & McKee (2002) studied rotational effects on halting the protostellar 
accretion onto metal-free stars using a simplified treatment. They point out the difficulty 
of forming very massive objects (M^^lOO Msun) because of the flow reversal in the polar 
directions. However, we believe that this conclusion is uncertain since accretion can still 
proceed through the circumstellar disk. Thus, the problem of making very massive stars is 
still open. 

6. Conclusions 

We have studied the evolution of extremely low metaUicity {Z/Zq = — 5 x 10~^) 
protostars accreting at very fast rates Mace ~ 10~^ — 10~^ Mq yr~^, as is expected to occur 
in the first star forming sites of cosmological structures. The main results can be summarized 
as follows: 

• The earhest stages of the protostellar growth are quahtatively the same independent 
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of the details of the evolution of the mass accretion rate, or of the presence of trace 
abundances of heavy elements. These phases correspond to the adiabatic accretion, 
the propagation of the luminosity wave, and the Kelvin-Helmholtz (KH) contraction. 
However, the mass range covered in each phase does depend on Mace the 
parameter space explored here extends up to ~ 60 Mq. 

• For more massive protostars, the evolution follows two fundamentally different branches, 
depending on the instantaneous value of the mass accretion rate. If Mace < Mc^t ^ 
4 X 10~^Mq yr~^, protostars relax to ordinary H-burning objects while still accreting. 
In the opposite case. Mace > Merit, the protostellar luminosity reaches the Eddington 
limit before the onset of nuclear burning at the center. This phase is soon followed by 
a rapid expansion of the radius, possibly leading to reversal of the accretion flow. 

• The two regimes yield upper mass limits that differ considerably: for Mace ^ Merit, 
the maximum mass is M^ax ^ 300 Mq and decreases with Mace; in the opposite case, 
Mmax ^ 300 Mq, allowing the formation of very massive objects. 

• Under a more realistic case of time dependent accretion, such as that proposed by Abel 
et al. (2002) where an initially high Maee begins to decrease only at high protostellar 
masses (> 60 Mq), the evolution resembles the case Maee < Merit. Therefore, the 
star is expected to grow in mass until the accretion rate drops substantially and the 
protostar joins the main sequence. At this time, the expansion of an ionized region 
would likely stop further accretion (> 460 Mq; Omukai & Inutsuka 2002). This limit 
is similar to that imposed by the amount of mass that can be accreted by the star 
during the entire main-sequence lifetime (600 Mq; Abel et al. 2002). 

• In all cases considered, deuterium burning and the pp-chain do not produce enough 
energy to counteract the effect of Kelvin-Helmoltz contraction. It is only with the 
beginning of the CN-cycle that nuclear burning has significant consequences for the 
protostellar evolution. 

• A slight amount of metals {Z/Zq ~ 10~^ — 10"^) favors continuing accretion by in- 
creasing the value of the critical accretion rate below which accretion can continue 
unimpeded. Thus, a non-zero metallicity favors the formation of very massive proto- 
stars. Above a threshold metallicity of Z ~ 0.01 Zq, the situation is reversed because 
of the feedback of radiation pressure on dust grains, as in the standard solar metallicity 
case. 

• Almost throughout the protostellar phase, the accreting envelope is optically thick. 
Although the stellar surface is as hot as ~ 10^ K, the photospheric temperature remains 
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limited to ~6000 K. Only after the protostellar mass exceeds > IOOM0 and the density 
in the envelope decreases substantially (~ 5 x 10~^^g cm~^ at the photosphere), the 
photospheric temperature begins to increase. In the case of the Abel et al. evolution, 
the photosphere coincides with the stellar surface for protostellar masses above SOOMq. 
Beyond this value, the hot stellar surface becomes optically visible. 

• Because of the high effective temperatures (~ lO^K), population III stars have been 
considered an efficient source of H and He ionizing photons, particularly in the case 
of star clusters consisting only of primordial ZAMS stars (Tumlinson & Shull 2000; 
Bromm ct al. 2001b; Schaercr 2002). However, during the main accretion phase, the 
ionizing photons arc degraded to optical photons in the precursor and cannot create a 
standard HII region. If the accretion phase is prolonged, the total emissivity of ionizing 
photons might be reduced by a non-negligible fraction with important consequences 
on the history of the re-ionization of the universe. 
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Figure Caption 
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Fig. 1. — Mass-radius relations for metal- free protostars evolving with an accretion rate 
-^acc — 1/4, 1/2, 1, 2Mfi(j (from top to bottom). The dashed line is for Mace — -^fid, where 
the fiducial value is Mgd = 4.4 x 10~^Mq yr"-*^. The filled circles indicate the onset of 
H-burning via the CN-cycle. 
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Fig. 2. — The internal structure of a protostar accreting at Mgcc = l/4Mfid = 1.1 x 
lO^'^M0/yr as a function of mass. Top panel: Evolution of the nuclear burning region 
and of the convective zone as a function of the relative mass M/M^. Regions where the 
energy generation by D and H burning exceed 10% of the average energy generation rate 
L^/M^ are shown as hatched areas. The extent of the convective zone is displayed by the 
dotted area. The sohd curve labeled Pr = Pg is the locus where the radiation pressure equals 
the gas pressure. Radiation pressure dominates to the right of this curve. The location of 
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the temperature peak is also indicated by the dotted line labeled Tmax- Bottom panel: The 
evolution of temperature and density at the center (solid lines) and at the position of the 
peak temperature (dotted lines). The CN abundance in the convective core is also shown. 
Each variable is normalized to 10^ K, 1 g cm~^, and 10~^, respectively. 
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Fig. 3.— Same as Fig. 2, but for M^cc = ^fid = 4.4 x IQ-^Mo/yr. 
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Fig. 4. — The evolution of the interior luminosity (solid) and the contribution due to deu- 
terium burning (dotted), the p-p chain (short-dashed), and the CN-cycle (long-dashed) for 
Mace = l/4Mfid = 1-1 X IQ-^Mq/yy {top panel) and for M^^ = Mgd = 4.4 x IQ-^Mq/yt 
{bottom panel). 



-27- 




Fig. 5. — Evolution of the interior luminosity-to-mass ratio as a function of mass for different 
values of the ratio Macc/Affid- The dashed horizontal line corresponds to the Eddington limit 
due to electron scattering. 
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Fig. 6. — Evolution of the core radius (solid) and the Eddington radius (dashed) for values 
of the accretion rate lower {top panel) and higher {bottom panel) than the fiducial one. 
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Fig. 7. — The evolution of the photospheric radius for two cases with l/4Mfi(j and Mfid, 
respectively. The dashed lines represent the evolution of the core radius. 
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Fig. 8. — HR diagram for primordial protostars. For comparison, we also show the locus 
of the metal-free ZAMS stars as computed by Marigo et al. (2001) for < lOOM© and 
by Bromm et al. (2001b) for higher masses. The dotted lines show the temperature and 
luminosity at the stellar surface. 
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Fig. 9. — The evolution of the protostellar radius with the time dependent accretion rate of 
Abel et al. (2002). The filled circle marks the onset of H-burning. The thin dashed lines 
represent the evolution for constant Mace = 1/4, 1/2, 1, and 2Mcrit- 
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Fig. 10. — Run of the photospheric radius in the time dependent accretion rate models. The 
protostellar radius is shown by the dashed hne. 
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Fig. 11. — The protostellar mass-radius relations at fixed Mace = l/2Mfid for different 
mctallicities Z = 0, lQ-^{= 5 x IQ-^Zq), 10-^(= 5 x IQ-^Zq). The filled circles mark the 
beginning of the CN-cycle. 
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Fig. 12. — The protostellar mass-radius relations computed for Mace = -^fid = 4.4 x 10 ^Mq 
yr"^ at different metallicities. 



